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This paper concerns itself with the analysis of a type of periodically 
switched network that might be used in time multiplex systems. The econom- 
ics of the situation require that the ratio of the switch closure time t to the 
switching period T be siyiall. ['sing this assumption, the analysis is performed 
by successive approximations. More precisely the zeroth approximation to 
the transmission is obtained from a block diagram analogous to those used 
in sampled servomcchanisms. From the convergence proof of the successive 
approximation scheme, it follows that when t/T is small, the zeroth approxi- 
mation is very close to the exact transmission. A discussion of some examples 
is included. 

CONTENTS P, lK c 

I. Introduction 1 404 

II. Description of the system 1404 

III. Method of solution 1407 

IV. The zeroth approximation 1408 

4.1 Introduction 1408 

4.2 Description of the block diagram 1410 

4.3 Analysis of the block diagram 1411 

V. Transmission loss 1413 

VI. A simple example 1414 

VII. Numerical examples 1416 

VIII. The successive approximation scheme 1417 

8.1 Preliminary steps 1418 

8.2 Matrix description of the successive approximations 1420 

5.3 Convergence proof 1421 

IX. Modification of the block diagram to improve the zeroth approxima- 
tion . . 1422 

X. Conclusions 1423 

Appendices 

I. Analysis of the resonant circuit 1423 

II. Study of the limiting case '/'->() 1425 

III. Zeroth approximation in the case where N\ is not identical to Nj 1425 

IV. Derivation of equation (24) 1426 

1403 



1404 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1957 

I. INTRODUCTION 

One main contributor to the cost of transmission circuits is the trans- 
mission medium itself. Thus it is important to share the transmission 
medium among as many messages as possible. One possible method is 
the frequency multiplex where each message utilizes a different frequency 
band of the whole band available in the medium. An alternate method 
is the time multiplex where each message is assigned a time slot of dura- 
tion t and has access to that time slot once every T seconds. It is obvious 
that the economics of the situation requires that r be as small as possible 
and T as large as possible so that the largest possible number of messages 
are transmitted over the medium. For this very reason the analysis of 
periodically switched networks is of special interest in the case where 
t/T is small. 

W. R. Bennett 4 has published an exact analysis of this problem without 
any restrictions either on the network or on the ratio t/T. It is believed, 
however, that the analysis presented in this paper will, in most practical 
cases, give the desired answer with a considerable reduction in the 
amount of calculations. The simplification of the analysis is mainly a 
result of the assumption that t/T is small. 

First the successive approximation method of solution will be discussed 
in general terms. Next it will be shown that the zeroth approximation 
to the transmission through the network can be obtained from the gain 
of a block diagram analogous to those used in the analysis of sampled 
servomechanisms. The nature of the zeroth approximation is further 
clarified by some general discussion and some examples. Next it is shown 
that the successive approximations converge. The convergence proof then 
suggests some slight modifications of the block diagram to obtain a more 
accurate solution. 

II. DESCRIPTION OF THE SYSTEM 

The system under consideration is shown on Fig. 1. It consists of two 
reactive networks Ni and N 2 connected through a switch S which is it- 
self in series with an inductance l. Ni is driven at its terminal pair (1) 




Fig. 1 — System under consideration. 
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Fig. 2 — Resonant circuit. 



by a current source h which is shunted by a one ohm resistor. JV 2 is also 
terminated at its terminal pair (1) by a one ohm resistor R L which is the 
load resistor of the system. The switch S is periodically closed for a dura- 
tion t. The switching period is T. Thus if the switch is closed during the 
interval (0, t) it will be closed during the intervals (nT, nT + t) for 
n = 1, 2, 3, • • • . The inductance I is selected so that the series circuit 
shown on Fig. 2 has a resonant frequency f r = 1/2t; i.e., the time r 
during which the switch is closed is exactly one-half period of the circuit 
of Fig. 2. 

The switch & acts as a sampler and, as a result of the well-known modu- 
lating properties of sampled systems, the sampling period T must be 
chosen such that the frequency 1/2 T is larger than any of the frequencies 
present in the signals generated by h . Furthermore, in order to eliminate 
all the sidebands generated by the switching, A r 2 must have a high in- 
sertion loss for all frequencies above 1/2T cps. 

In the analysis that follows networks Ni and N 2 will be assumed to be 
identical: it should, however, be stressed that this assumption is not 
necessary for the proposed method of anafysis.* This assumption is 
made because in the practical situation which motivated this analysis 
iVi and iV 2 were identical since transmission in both directions was re- 
quired. 

In order for the system under consideration to achieve the maximum 
degree of multiplexing, the closure time r of the switch will be taken as 
small as practically possible and the switching period T as large as pos- 
sible (consistent with the bandwidth of the signals to be transmitted). 
As a result the ratio t/T is very small, of the order of 10~ 2 or less in prac- 
tical cases. Consequently the resonant frequency f r of the series resonant 
circuit shown on Fig. 2, is many times larger than any of the natural 
frequencies of JV"i and iV 2 • 



* The modifications required for the case where iVi is not identical to JV« are 
given in Appendix IV. 
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The problem is to determine the relation between Et , the voltage 
across R L , and h . 



III. METHOD OF SOLUTION 

Let us first write the equations of the system. Obviously the equations 
will depend on the exact configuration of the networks N x and N 2 . For 
simplicity we shall write them for the case where Ni and ATj are dissipa- 
tionless low-pass ladder networks. As will become apparent later this 
assumption is not essential to the argument. What is essential, however, 
is the fact that both Ni and JV 2 should start (looking in from the switch) 
with a shunt capacitor C and a series inductance L„ , the element value 
of L„ being much larger than C. Using a method of analysis advocated 
by T. R. Bashkow, 5 we obtain, for the network of Fig. 3, the equations: 

L,^= -RU - v, + ftio 



& Tt = ii - H 



h 



n dv n _ . 

Un —T— — ln_ 1 1 n 

L n ^ = v n -e 2 (l.a) 

dt 

C 1 ^; = u-i r A(t) (l.b) 
at 

C^j = [* - e 3 ]A{t) (l.c)U 
at 

c^ = i,A(t) -i,: d.d) 



T di n ' I 

Ln-7- = e 3 - v n 
dt 



n dih . i ■ i 
at 

T di\ l D • ' 

dt 



> h 



(0 
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where 

Mt) = £ W(t - kT) - u(t - kT - r)], (2) 

fc— O0 

with u(t) = 1 for t > 0, and u(t) = for t < 0. 

This system of linear time varying equations may be broken up irto 
three sub-systems I a , R and h . It is this subdivision that suggests a 
successive approximation scheme that will be shown to converge to the 
exact solution. 

The zeroth approximation is obtained as follows: when the switch is 
closed, i.e., A(l) = 1, the resonant current i T is much larger than the cur- 
rents i n and i,' . Thus, during the switch losure time, i n and i n ' are neg- 
lected with respect to i r in (Lb) and (l.d). Hence when A(t) = 1 the 
system R may be solved for i r (t), e 2 (t) and e s (l) in terms of the initial 
conditions. The resulting function e 2 (i) and given function i (t) are then 
the forcing functions of the system /„. The other function e 3 (t) is the 
forcing function of the system h. Under these assumptions, the periodic 
steady-state solution corresponding to an applied current i Q (t) = Toe'"" 
is easily obtained. 

The zeroth approximation will be distinguished by a subscript "0". 
Thus iro(t) is the (steady state) zeroth approximation to the exact solu- 
tion i r (t). 

The first approximation will be the solution of the system (1), pro- 
vided that during the switch closure time the functions i„(t) and i n '{t) 
in (l.b) and (l.d) are respectively replaced by the known functions 
i n0 (t) &ndi„ '(t). And, more generally, the (k + l)th approximation will 
be the solution of (1) provided that during the switch closure time, the 
functions i n (t) and i n '(l), in (l.b) and (l.d), are respectively replaced by 
the known solutions for i n (t), and i n '(t) given by the fcth approximation 
It will be shown later that this successive approximation scheme con- 
verges. Let us first describe a simple method for obtaining the zeroth 
approximation. 

IV. THE ZEROTH APPROXIMATION 

4.1 Introduction 

The problem is to obtain the steady-state solution of (1) under the 
excitation i (t) = I e ut . Using the approximations indicated above, 
during the switch closure time (that is when &(t) = 1) the system R 
becomes 
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C Tt = - ir{t)m > (3) 

t d £ = \e, - e 3 }A((), (4) 

at 

Cj* = i r (t)A(t). (5) 

Differentiating the middle equation and eliminating dc 2 /dt and de 3 /dt 
we get f or ^ t < r : 

W = ~ Jc iMt) + \ [e2il) ~ e3{mo (0) 

in which we used the notation v'{t) for the Dirac function and the knowl- 
edge that 

^ . 5(0 _ m _ r) (7) 

Eciuation ((>) represents the behavior of the resonant circuit of Fig. 2 
for the following initial conditions: 

*V(0+) = 0, (8) 

tf*V(0+) c 2 (0) - e 3 (0) 



dt I 



(9) 



In Appendix I it is shown that the resulting current i r (t) is, for the in- 
terval ^ I < T, 

ir(t) = C[c 2 (0) - e 3 (0)] Sl (0, (10) 

where 

— - sin — = - too sin u) t for ^ t < r 
s,()=| 2t - 2 (11) 



, elsewhere 
with 

7T . /2 



■;"V<e- (12) 

Thus the zeroth approximation to the exact i r (t) is given for the interval 
^ t ^ T by 

M(0 = C[e 2 (0) - c 3 (0)]si(0. (13) 
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We shall now show that the zeroth approximation may be conveniently 
obtained from the block diagram of Fig. 4. 

4.2 Description of the Block Diagram 

All the blocks of the block diagram are unilateral and their correspond- 
ing transfer functions are defined in the following. Capital symbols repre- 
sent £-transform of the corresponding time functions, thus 7 (p) is the 
43-transform of i (t). 

Referring to Fig. 1, 



ZM = 



E,(p) 

Up) 



/ r =0 



Thus Z i2 (p) represents the transfer impedance of Ni when its output is 
open-circuited (i.e., l T = 0). Since Ni and N 2 are identical we also have, 
from R L = 1 and reciprocity, Z u (p) = Ei/I r , where I r is the cur- 
rent entering N? . 

The impulse modulator is periodically operated every T seconds, 
and has the property that if its input is a continuous function f(t) its 
output is a sequence of impulses: 

Zf(t)8(t-kT). 



The transfer function Si(p) is defined by 



Si(p) = jbMOI = -a 



Wn 



cosh ~ e p . 



;i4) 



p~ + uo 



Let Z(p) be the driving point impedance at the terminal pair (2) of Ni . 
It is also that of N 2 since Ni and N 2 are assumed to be identical. 

Let V(p) be the output of the first block, then, by definition, V(p) = 
Z n (p)h ■ Let v(t) be the corresponding time function. The voltage v(t) 
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Fig. 4 — Zeroth-approximation block diagram. 
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is the output voltage of N x , when Ni is excited by the current source 
/o and the switch S remains open at all times. 

4.3 Analysis of the Block Diagram 

For simplicity, suppose that the system starts from a relaxed condi- 
tion (i.e., no energy stored) at / = 0. Let z{t) = £~ [Z(p)]. Considering 
the network Ni as driven by i and i T o , it follows that the voltage e->(t) 
shown on Fig. 3 is given by 

cM = v(t) - [ iM)z(t - t') dt'. (15) 

Similarly 

e 3 o(0 = I Ut')z(t -t')dt'. (16) 

Thus 

eM - c, a (t) = v(t) - 2 [ iM)z(l - t') dt'. (17) 

Jo 

These equations have been derived by considering Fig. 1. They could 
have been also derived from the block diagram of Fig. 4 as follows: let 
/r(»(p) be the output of G\S*i(p). As a result, the output of the block 
2Z(p) is 2Z(p)I r o(p). When this latter quantity is subtracted from V(p) 
one gets V(p) — 2Z(p)/ rl) (p), which is the ^-transform of the right-hand 
side of (17). Referring to the block diagram it is also seen that this 
quantity is the input to the impulse modulator. 

Thus we see that if 7h>(p) is the output of CSi(p), then the input of 
the impulse modulator is e 2 o(0 — e 3 o(t) by virtue of (17). If this is the 
case the output of CSi(p) is given by C[e 20 (0) - e 3 o(0)]si(0, for g t < T, 
which, according to (9), is i r o{t). 

Thus the block diagram of Fig. 4 is a convenient way of obtaining the 
zeroth approximation to the periodic steady-state solution. 

In order to use the techniques developed for sampled data systems, • 
we introduce the following notation. 2 If f(t) = ^^[Fip)], then we define 
F*(p) by the relation 

F*(P) =4 £ /<Xp+>co s ), (18) 

1 H — 00 

where 

co. = ^ . (19) 



1412 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1957 

If /(0+) is denned by lim/(e 2 ), then, provided /(0+) = 0,* 



F*(p) = £ 



t, fit) S(t - nT)]. 

n=0 J 



(20) 



Going back to the system of Fig. 4 we get" 



L40(P) ~ 1 + 2C[Si(p)Z(p)j* ' U ^ 



and 



Up) = 1 + 2QSMZW (22) 

where according to the notation defined by (18) 

[Z„(p)/o(p)]* = ^ Z ZAp + jnuMp + ./no,,), 

1 H=— 00 

1 + " 

[Si(p)Z(p)]* = 1 £ 5i(p + jm*)Z(p + ./no,.). 

i n=— oo 

It should be stressed that (21) and (22) are not valid when t is made 
identical to zero. When t = 0, Si(p) = 1 for all p's and since Z{p) ~ 
1/Cp as p — * «> the time function whose transform is Z(p)Si(p) is differ- 
ent from zero at I = 0. In such a case (20) does not hold. From a physi- 
cal point of view, the feedback loop of Fig. 4 is unstable when t is identi- 
cally zero since an impulse generated by the impulse modulator 
produces instantaneously a step at the input of the impulse modulator. 
This step causes an instantaneous jump in the measure of the impulse 
at the output of the impulse modulator and so on. In short the feedback 
loop is unstable. 

It should be pointed out that if the power density spectrum of h is 
zero for frequencies higher than o>,/2, (21) reduces to 

S,l C&(P)*»'« for |p|<-=. (23) 

h n + 2C[S l (p)Z(p)}* ' p ' 2 

For certain applications it is convenient to rewrite (21) in a slightly 

* When/(0), as defined above, is different from zero, (20) should be replaced 
by 

£ |~£ fit) bit - nTU = F*(p) + \ /(0+). 
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different form. Advancing the time function si(t) by r/2 seconds, one 
gets the function s () (/) which is even in t. As a result its transform So(p) 
is purely real, that is, 

2 

So(p) = , a> ° ■ cosh^. 

P~ + O) ^ 

From an analysis carried out in detail in Appendix IV we finally obtain 

F , x _ [Z 18 (p)/o(p)]* < S,(p)Z 18 (p) r . 

2[Z(p)So{p)r 

It should be pointed out that (23) is still valid when r = 0. Equations 
(20) and (23) give the zeroth approximation to the gain of the system 
for any driving current io(t). 

In many cases it is sufficient to know only the steady-state response 
7?.io(p) to an input i (t) = I e"* ' . The response E i(s {p), as given by (23) 
[or (20)] includes both transient and steady -state terms. Since h(p) = 

V- equation (24) gives 

P - jo>o 

(l £ Z 12 (p + inoO . /o J) S (p)Z 12 (p) 

*•<«-- 2[Z(p)&(p)]* "• (2o) 

Since neither S n (p) nor Zu(p) have poles on the imaginary axis, the 
steady state includes only the terms corresponding to the imaginary 
axis poles of the summation terms. Thus the steady-state response is of 
the form 

g An e /(»„-»..)« f 



where, from (25), 

2 £ Zljwo + j(* - nMSoljm + j(* - «)«.] 



, _ IoZi2(jcoo)So(jua — jnu s )Zn{jm — jnoi s ) 

- ln +5 ■ . (26) 



k^x 



V. TRANSMISSION LOSS 



A practically important question is to find out a priori whether a 
switched filter necessarily introduces some transmission loss. 

The following considerations apply exclusively to the zeroth order 
approximation. It will be shown that assuming ideal elements, the trans- 
mission at dc may have as small a loss as desired. 
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By transmission at dc we mean the ratio of the dc component of the 
steady state output voltage to the intensity of the applied direct current. 
Thus we refer to (26) and set co„ = and n = 0. Suppose the lossless 
networks Ni and N» are designed so that their transfer impedance Z n 
is of the Butterworth type, that is 

I MM I 2 = tz^> 

where for our purposes M is a large integer. 

In the following sum, which is the denominator of (20) when am = 
n = 

2 £ Z(jka,)S (jku 8 ), 

k=-x 

(where u a > 2 since the cutoff of the networks N occurs at w = 1), the 
terms corresponding to values of k ^ will make a contribution that 
vanishes asM-* oo . This is a consequence of the following facts: 

(a) Re[Z(jJfe»,)] = | ZiaO'fao.) 1 2 , since the networks tf x and JV 2 are 
dissipationless. Hence for A; ^ and as M — > <*> Re[Z(jko),)] — > 0, 

(b) hn[Z(jko> a )] = -Im[Z(-jku,)], 

(c) S {jo>) is real. 

Thus the imaginary part of the products Z{jku a )So(jkta,) cancel out and 
the real part (for fc ^ 0) decreases exponentially to zero as M — » » . 
Hence for sufficiently large M the denominator of (20) may be made as 
close to two as desired. 

It is easy to check that the numerator of (20) reduces to h , the in- 
tensity of the applied direct current. Therefore the ratio of A , the dc 
component of the output voltage to h may be made as close to one-half 
as desired. 

VI. A SIMPLE EXAMPLE 

Since the approximate formulae derived in Section IV are somewhat 
unfamiliar it seems proper to consider in a rather detailed manner a 
simple example.* 

Consider the system of Fig. 5. Assume that the current source applies 
a constant current to the system and assume that the steady state is 
reached. For simplicity let hR = E. 

The steady-state behavior of the voltages e 2 (t) and c 3 (() = c 4 (t) is 

* In addition, the limiting case of the sampling rate — > », i.e., T — > 0, is treated 
in Appendix II. 
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1 



~ In R<" C: 



:C ->R e. 



Fig. 5 — A simple example. 



V,+A r 




Fig. 6 — Waveforms of the network of Fig. 5. 

shown on Fig. 6. It is further assumed that the duration r during which 
the switch is closed is negligible compared to T, the interval between 
two successive closures. 

Let e» be the average value of the steady-state voltage e\{t). Thus £4 
is equal to A as given by (26) with w = n = 0, namely, 

Z 12 2 (o) So(Q) 

64 — -*0 +=o 

2 J] Z(jku>,) SoO'A-co.) 



In this particular case 

Z(p) = Zn(p) = 



R 



C~ l 



1 + pRC i • 



Since we assume r to be infinitesimal S (p) and S\(p) may be con- 
sidered equal to unity over the band of interest. Using the expansion 
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1 °° 

coth Z = - + /2 



2z 



z i z 2 +nV 
we obtain 

i +■» c~ l 1 

Z ^4g , , ■ , . 1 = 2C C ° th 
(p + jnojs) + ^ 



P + ^)I- 



Hence 



Thus finally 



% 2 ZOAtoO = 2Z*(p) 



P=0 = ^ coth Gy • 



/o^ 2 „ tf C 



c = -, , \ = E-— / t \ (27) 

ricoth(^) coth '2«cJ 



This last result obtained from the theory developed above is now going 
to be checked directly. Referring to Fig. 6, where the notation is denned, 
and noting the periodicity of the boundary conditions, we get 

(7x + A)e- {TIRC) = V, , 

(JE - V 1 )e- {TIRC) = E - (Vi + A). 
Noting that e 4 (t) = (Vi + A)e~ t/BC , and solving for Vi and A we fin- 
ally get 

-tIRC 
C4«) - ' 



1 + e- T ' RC ' 
By definition 

_ i r , A .. erci- r T,RC 



or 



1 f rA ,. ERC l-S "~ 



RC 1 

ei — E m / T \' 

coth WJ 

This last equation checks with (27). 

VII. NUMERICAL EXAMPLES 

Consider the network of Fig. 7. The cutoff of both iVi and N 2 occurs 
at co = 1. In view of the sampling theorem good transmission requires 
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2.0 




Fig. 7 — Computed transmission loss. 

that the signal be sampled at a rate at least twice as large as its highest 
frequency component. Since the cutoff occurs at co = 1, the sampling 
angular frequency should at least be equal to 2. For illustration purposes 
we have taken g>„ = 2. 07 and w„ = 5 for the angular sampling frequency 
The value w» = 2.67 corresponds to a cutoff at 3 kc and a sampling rate, 
of 8 kc. The ratio t/T was taken to be 1/125. The transmission through 
the switched network as given by the zeroth approximation is shown for 
both cases on Fig. 7. 

As expected the transmittance of the switched filter gets closer to that 
of an ordinary filter as the switching frequency increases. 



VIII. THE SUCCESSIVE APPROXIMATION SCHEME 

The ideas involved in the successive approximation scheme are simple 
and straightforward. One point remains to be settled, namely the con- 
vergence of the procedure. 

We shall assign a subscript 1 to the correction to be applied to the 
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zeroth approximation in order to obtain the first approximation. Thus 
adding i rl (t) to *Vo(0 we get the first approximation ^(0 + i T i(t). More 
generally the A-th approximation is ^C—ntruft). The procedure will con- 
verge if, in particular, the infinite series sX-t(im(t) converges. 

8.1 Preliminary Steps 

(a) Let us normalize the frequency (and consequently the time) so 
that the switching period T is unity. Since the networks JVi and N* must 
have high insertion loss for to > %(2ir/T) = t, the pass band of Ni and 
JV 2 must be the order of 1 radian/sec. As a result the element values of 
the capacitor C and the inductance L„ (see Fig. 3) are also 0(1). 

(b) For the excitation i = e ,al , the zeroth approximation derived 
above may be written in terms of Fourier components: 

+» 

troW = e 2-, 1 ro,ke , 

k=-oo 

+00 

= e 2-> in0,k6 

fc=— oo 

Let ho denote the complex conjugate of i r o , then 

+ 00 

iro(t)lM = 2 Ir0,kIrO.ie 2Tl{k ~ tU . 
k,(=-ao 

Since the functions e 2 ' iki [k = • • • — 1, 0, 1, • • • ] are orthonormal over 
the interval (0, 1) and form a complete set, 3 we have from Bessel's equal- 

ity: 3 

.1 +00 

i ra (t) \'dt = £ I / r0lt | 2 = N(Iro), 
)0 fc=_oo 

where N(I r o) denotes the norm of the vector I r0 which is defined by its 
components I, ,k(k = •'• — 1| 0. +1 ■ * * )• Similarly, 

' | ino(0 | 2 dt = S I /-.* I 1 - N(T n0 ). 



f 

Jo 



J 



(c) Since the switch is periodically closed we shall be interested in 
the Fourier series expansion of A(/): 

A(0 =u(t) -n(t-r) =EV K ", 

where A = r and 

Z I A fc | 2 = r. 

— 00 
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Since t/T « 1, and since the frequency has been normalized so that 
T = 1, we have t« 1. 
Using the convolution in the frequency domain, we have 



U0A(0 = Z ( Z A 

A - =— oo \o=— oo 



:-Jna, a I C 



f(u+2»fc)( 



If we introduce the infinite matrix G defined by 

G ik = A,-_ fc (t, A; = - =o, • • • , -1, 0, +1, ■ • • oo), 

the convolution may be represented by the product, GI n o , where /„o is 
the vector whose components are 7„ , a(A' = • • • —1, 0, 1, • • •)• 

(d) Considering the network shown on Fig. 3, let E{p) be the ratio of 
I n '(.p) to I r (p). Taking into account the assumed identity between ATi 
and N 2 it follows that 



+/.(p) 



I rip) 



xS - *» 



Using the system of (1) and, for example, by Neumann series expan- 
sion of the inverse matrix, we get 

E(v) = — — + • • • 

(e) Considering now the effect of i„(t) and i n '(t) on i r (t), (42) of 
Appendix I gives I r (p) as a function of I n (p) and I n '(p). In the present 
discussion where we are interested in the steady state of i r (t) it is es- 
sential to keep in mind that since the switch opens at t = r, the memory 
of the resonant circuit extends only over an interval < t ^ r. To take 
this into account we must modify the factor (w 2 /2)/(p 2 + w 2 ) of (40), 
because the impulse response (which represents this memory) must be 
identically zero for t > t. The resulting new expression is 

2 

F(p) - -r^—i e-" T ' 2 \e p,n + e- pT ' 2 ] , 
p- + wo- 
or 

F(p) =_^L_ e -^ 2 C o S h^ T . 

p- + coo 2 

Since the time function whose transform is F(p) is non-negative for all 
t's and since /<'(0) = 1 , it follows that 

\F(Ju) | ^ 1. (28) 
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8.2 Matrix description of the successive approximations 

From the developments of Section IV, we know i, (t), i n o(t) and i n o'(t) 
or what is equivalent, the vectors 7 r0 , I n o and /„</. The first approxi- 
mation takes into account the effect of i n0 (t) and i, M '{t) on i T (t). [See 
equation (l.c) and (l.d)]. The time functions i n0 (0 and i„o'(0 affect the 
system R only during the interval (0, t). Therefore we must consider the 
vector G{I„o + /no') which corresponds to the excitation of the resonant 
circuit. Since the opening of the switch after a closure time r forcibly 
brings i r (t) to zero we have 

In = GfG(/„o + U (29) 

where the matrix G has been defined above and the matrix F is a diago- 
nal matrix whose diagonal elements F k (k = • • • — 1, 0, +1 • ■ •) are 
defined by F k = F(ju>, + j2irk). Note that (28) implies that | F k \ £ 1 
for all fc's. It should be kept in mind that I r0 + In is the first approxima- 
tion to the exact I r (p)- 

The next iteration is obtained by first taking into account the effect 
of In on the rest of the network: 

/nl = E In ' (80) 

Inl' = E In , 

where E is a diagonal matrix whose elements E k (k = ■ ■ • , —1, 0, +1, 
■■■) are defined by E k = E(ju° + 2irfcj)j and then the effects of I n2 
and I, ,2 on 7, , that is, 

/,, = GFG{I ni + /„/), (31) 

combining (30) and (31), / r2 = 2GFG E In . A repetition of the same 
procedure would lead to I r3 = 2GFG E Li , and in general I rn + 1 = 
2GFGEI rn . 

Since the nth approximation to I r (p) is given by the sum X^U/rt , 
the successive approximation scheme will be convergent only if the series 

co 



Irk 

A-=0 



converges. This will be the case if and only if the series 

[1 +2GFGE + ■■■ + {2GFGE) n + ■ - -}In (32) 

converges. 
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8.3 Convergence Proof 

Consider a vector X of bounded norm corresponding to a time func- 
tion x(t) having the property that x(t) = for t ^ t ^ T and x(t) 9^ 
for < t < t. In the above scheme, the vector X would be l rn . Let us 
define the vectors Y, Z, U and V by the relations 

Y = EX, (33) 
Z = GY, (34) 
U = FZ, (35) 

V = 2GU, (36) 
hence 

V = 2GFGEX. (37) 

We wish to show that N(V) ^ aN(X) with a < 1, since these inequali- 
ties imply that the infinite series (32) converges. 

Since (a) Ni and N 2 are low-pass filters with cutoff ^ ir radians/sec, 
(b) E(p) = 1 for p = 0, (c) E(p) « \/L n Cp for p » 1, only a few of 
the .&ys will be of the order of unity In most cases E-\ , E , Ei will be 
smaller than unity, thus, 

N(Y) ^ N(X). (38) 

In view of the pulsating character of x(t) the power spectrum of x(t) 
is almost constant up to frequencies of the order of t/t radians/sec. 
Because of the low-pass characteristic of E(p), the function y(t) associ- 
ated with the vector Y is smooth in comparison to x(t), thus from (34), 

N{Z) = [ I 2(0 | 2 dt = f I y(t) p dt = arN(Y), 

JO JO 

where a = 0(1). 

Since | F k | £ 1 for all fe's, from (33), N(U) ^ N(Z), hence N(U) = 
brN(Y) with b = 0(1). 

N(U) = brN(Y) with b = 0(1). 

From (36) we have 

N(V) = 2 f I u(t) | 2 <// g 2 [ I u(0 | 2 tf/ = 2N(U). 
Jo Jo 

Thus we finally get 

N(V) = 2b T N(Y) where b = 0(1), (39) 
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and since t « 1 we get from (38) and (39) N(V) = aN(X) with a < 1. 
Hence the convergence is established. 

IX. A MODIFICATION OF THE BLOCK DIAGRAM TO IMPROVE THE ZEROTH 
APPROXIMATION 

In principle it is possible to obtain a block diagram whose transmission 
characteristic is equal to the first approximation. In many cases it is 
not necessary to go that far. The first approximation takes into account 
the effect of the currents i n0 (l) and i nQ '(t) on the resonant circuit of Fig. 2. 
Since during the switch closure time the currents i„ a and i„o' cannot vary 
much, let us assume that they remain constant for the duration of the 
switch closure. 

Referring to the analysis of Appendix I and to (42) in particular, we 
see that the current i r is increased by 



8i n (p) = 



COo 



t.(0-) + in(0-) 



p 2 + coo 2 



2p 



or 



N (0 = C d ' (0 " ) = 4(Q " } (1 - cos cooO OIK, 
Defining S 2 (p) = ST^iO- - cos a Q t)[u{t) - u(t - t)]}, or 



Stip) = 



W|> 



(2p 2 + 



2 p(p 2 + coo 2 ) 2p(p 2 + coo ; 



coo J -. 




e4p) = 



c[z, 2 (p)I (p)]*S 1 (p)+c[pZ 1a (p)Io(p)]*S 2 (p) 
1+2C-f[z(p)S,(p)]*+ [pZ(p)S 2 (p)] 

Fig. 8 — Modified block diagram. 
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and recalling that the input of the impulse modulator of Fig. 4 is 
e 2 (0) — e : i(0), it becomes obvious that the modified block diagram should 
be that given by Fig. 8. The output of the modified block diagram is 
given by 2 

E , v = C[Z»(p)7o(p)]*gi(p) + C[pZ»(p)/ (p)]»&(p) 

W " i + 2C{[Z(p)S 1 (p)]* + [pZ(p)S 2 (p)}*} * 

X. CONCLUSION 

Let us compare the method of solution presented above with the more 
formal approach proposed by Bennett. The latter method leads to the 
exact steady-state transmission through a network containing periodi- 
cally operated switches. This method is perfectly general in that it does 
not require any assumption relative to the properties of the network 
nor to the ratio of t/T. As expected this generality implies a lot of de- 
tailed computations. In particular it requires, for each reactance of the 
network, the computation of the voltage across it due to any initial con- 
dition. The method presented in this paper is not so general because it 
assumes first that the ratio t/T is small; second the value of the induct- 
ance I is very much smaller than that of L„ (see Fig. 3). The result 
of these assumptions is that the system of time varying equations 
may be solved by successive approximations with the further advantage 
that the convergence proof guarantees that, for very small t/T, the 
zeroth approximation will be a close estimate of the exact solution. 

The zeroth approximation may conveniently be obtained by consider- 
ing a block-diagram analogous to those used in the analysis of sampled 
servomechanisms. Further the proposed method leads directly to some 
interesting results, for example, as far as the zeroth approximation is 
concerned, the dc transmission may be achieved with as small a loss as 
desired provided the lossless networks Ni and N 2 are suitably designed. 
Another advantage of the proposed method is that the simplicity of the 
analysis permits the designer to investigate at a small cost a large num- 
ber of possible designs. 

Finally it should be pointed out that this approach to the solution of 
a system of time-varying linear differential equations may find applica- 
tions in many other physical problems. 

Appendix I 

ANALYSIS OF THE RESONANT CIRCUIT 

Consider the resonant circuit of Fig. 2. Suppose that at t = 0, the left- 
hand capacitor has a potential e 2 (0) and the right-hand capacitor has 
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a potential e 3 (0) and that at t = the current i r through the inductance 
I is zero. 

The network equation is 



'a 4 + H '"" : " ,a 



Now ir(0) = and diM/dt = [e 2 (0) - e 3 (0)]/£ Let 2/CC = to 2 , then 
di r (0)/di = u 2 C[e 2 (0) - e 3 (0)]/2. 
Using Laplace transforms, 



(p 2 + co 2 )I r (p) = piM + 



<fc(0) 
dt ' 



/r(p) - ^ [C2(0) " 63(0)] ' 



(40) 



p 2 + coo 2 ' 



hence 



and 



■(A ^ g 2 (0) - e 3 (0) • , 

lr(i; = O) C jr sin CO t 



(41) 



qit) = f\ {tUt = c e M^M 



[1 — COS Wot]. 



If 2v/tao = 2t, i.e., r = vs/tC/2, which means that the duration of the 
switch closure is a half-period of the resonance of the tuned circuit, then 



. ( , t C[e,(0) - e,(0)] . irt 
U{t) = = sm T' 

T Z T 

gM _ CMP -^Wlfi -«,?' 



[l-cos^]. 



It is clear then that, during the period r, the charge transferred onto the 




Ir(P) = 



w 



In(p)+In(P) 



p' + a>< 



Fig. 9 — Resonant circuit excited by current sources /„ and /„' . 
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right-hand capacitor is q(r) = C[e 2 (0) — e 3 (0)] and as a result at time 
t = t the right-hand capacitor has a voltage e 2 (0) and the left-hand 
capacitor has a voltage e 3 (0). Considering now the network of Fig. 9, 
the equation is 

*% + §+ =^^«) +i/(01. 

Assuming all initial conditions* to be zero we get, 

'>(p) = -rT-> Up)+ 2 ! "' (P) - (42) 

p L + oi l 2 

Appendix II 

STUDY OF THE LIMITING CASE T — * 

We expect that if the sampling period T — ► 0, which is equivalent to 
stating that the sampling frequency co s — * °c , then the inductance I —> 
and as a result the voltage e 3 (0 will be infinitely close, at all times, to 
the voltage e 2 (t). Thus, in the limit, everything happens as if the termi- 
nal pairs (2) of iVi and N 2 were directly connected. In that case the gain 
of the system is 

Zn(p) 7 , x 

Wp) Zl2(p) ' 

as is easily seen by referring to the Thevenin equivalent circuit of Ni . 
Let us show that as T — > 0, (21) leads to the same result. First note 
that both Zuh and ZSi go to zero at least as fast as 1/p 2 for p — > co . 
Hence the summations in (21) reduce to the term corresponding to 
n = 0. Therefore, 

hiKv) - 1 + 2CIZ5J* Zl2 ? + 2CZS! ^ 2Z /o • 
Appendix III 

ZEROTH APPROXIMATION IN THE CASE WHERE Ni IS NOT IDENTICAL TO A^ 2 

Let, for k — 1,2; Ck be the shunt capacitor at the terminal pair 2 of 
Nk , Z k {p) be the driving point impedance of iV\. , and Zu (p) be the 
transfer impedance of Nk . In the present case the capacitors Ci and C 2 
are in series in the resonant circuit of Fig. 2. It can be shown that the 

* Their contribution has been found in (40). 
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charge exchanged during one-half period of the resonance is 

-^p^ MO) - «((»]. 

Ci -\- C2 
For the present case, (16) and (17) become 

i(t) = v(t) - I iro(r)zi(t - t) dr, 

J — cO 

hit) = J iroWteit - r) dr. 



C 2 l 



c,{ 



Following the same procedure as before we are finally led to the block 
diagram of Fig. 10 whose output is given by 



E A (p) = 



[Zrr(p)h(p)Y 



2C1C2 



SMZvTip) 



1 + 



2ClC2 |[Z 1 (p) + Z 2 (p)]5 1 (p)}* 



C, + Co 

Appendix IV 

THE DERIVATION OF EQUATION (24) 

Considering the method used in Section IV to derive the zeroth ap- 
proximation, it is clear that during the switch closure the voltages e 2 (t) 
and e 3 (i) vary sinusoidally, that is, 



e 2 (0 = e 2 (0) - 



e«(0) - e,(0) 



l~i * l ~] 

1 — cos — , 



eM . ,3(0)+ m^M [l_ COS * 



Io 





IMPULS 
MODULA1 


E 
OR 














^ 


2C,C2 S,(p) 


- — > 


z (2) 


+ s 




j 






J 


t 


















Zi 


(P)H 


-Z a (p) 





















Fig. 10 — Zeroth approximation: modified block diagram for the case where Ni 
and Nt are not identical. 
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Thus, it always happens that for t = t/2, i.e., at the middle of switch 
closure time, c 2 (/) — e 3 (t) = 0. 

Therefore if we consider the time function e 2 (t) — e 3 (t) we have for all 

fc' S (-oo, ... 0, ••• + °°), 

62 ( /CT + " 63 ( kT + 2~) = °- 

If, for simplicity of analysis, we assume that the switch is closed during 
the intervals -(t/2) + kT ^ t ^ +t/2 + kT, then for all k's, 

e 2 (kT) — e 3 (kT) = 0. 

Using (17), this condition implies that [V(p)]* - 2{I r0 {p)Z{p))* - 0. 

Now, remembering that iro(t) consists of a sequence of half sine waves 
whose shape is defined by s (t) (which is by definition identical to Si(t) 
except for an advance in time of t/2) it follows that I r o{p) = B(p)S (p), 
where B(p) is the £ -transform of the sequence of impulses whose measure 
is equal to the charge interchanged between Ni and iV 2 at each switch 
closure. Since [B(p)S (p)Z(p)]* = B(p)[S (p)Z(p)}*, then 

R( v [Zu(p)/o(p)l* 
W 2[S (p)Z(p)]*- 

From which it immediately follows that 

[Zi2(p)/o(p)]*£o(p) 



/ro(p) = 



and 



where 



E M (p) = 



2[S (p)Z(p)}* 

[Z»(p)/o(p)]»iSo(p)Zi a (p) 
2[S Q (p)Z(p)]* 



[S (p) Z(p)}* = ^ E £o(p + jn«.) Z(p + Jim.) . 

1 n = —oo 
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